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Abstract 

Hamiltonian lattice gauge models based on the assignment of the Heisenberg 
double of a Lie group to each link of the lattice are constructed in arbitrary space- 
time dimensions. It is shown that the corresponding generalization of the gauge- 
invariant Wilson line observables requires to attach to each vertex of the line a 
vertex operator which goes to the unity in the continuum limit. 


1 Introduction 

In a recent paper |lp I have considered Hamiltonian lattice Yang-Mills theory, based on the 
assignment of the Heisenberg double of a Lie group to each link of the lattice, in (1+1)- 
and (2+l)-dimensions. In the present paper having in mind possible applications of the 
construction to Chern-Simons models and gravity I discuss an analogous formulation of 
lattice gauge models on an arbitrary lattice or graph. In the case of the regular hyper- 
cubic space lattice the gauge models proposed are lattice-regularized Yang-Mills models 
in (d+l)-dimensions. 

As is well known there are two possible ways of lattice regularization of gauge theo¬ 
ries. In the approach of Wilson [2j] one considers Euclidean formulation of a model and 
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discretizes all space-time, thus replacing the model by some statistical mechanics model. 
In the Hamiltonian approach of Kogut and Susskind [3;j one considers a model in the 
Minkowskian space-time and introduces only space lattice remaining the time direction 
continuous. Then one places on each link of the lattice the cotangent bundle of a Lie 
group and on each vertex lattice Gauss-law constraints which are first-class constraints 
and generate gauge transformations and finally one finds a gauge-invariant lattice Hamil¬ 
tonian. Thus in the Hamiltonian formulation the Yang-Mills theory is replaced by some 
classical mechanics model with first-class constraints, the phase space of the model being 
the direct product of the cotangent bundle over all links of the lattice: Yl Hnks T*G. 

However one could ask oneself whether it is possible to place on each link another 
phase space and on each vertex other Gauss-law constraints and to get another lattice 
model which can be reduced to the continuous one in the continuum limit. In |l| I have 
shown that such a possibility does exist and is based on a phase space which is called 
the Heisenberg double D\ of a Lie group and is one of the basic objects in the theory 
of Poisson-Lie groups Si- Only (1+1)- and (2+l)-dimensional Yang-Mills models were 
considered in the paper. In the present paper a generalization of the consideration to the 
Yang-Mills theory in any space-time dimension is found. 

The plan of the paper is as follows. In the second section we remind some simple 
results from the theory of the Heisenberg double and introduce the notations used in 
the paper. In the third section we formulate a classical mechanics lattice gauge model 
with first-class constraints on arbitrary lattice or graph, the phase space of the model 
being the direct product of the Heisenberg double over all links: Ylu nks D^_. Then we 
show that in the case of the regular hyper-cubic space lattice in d-dimensions the model 
constructed is lattice-regularized (d+l)-dimensional Yang-Mills theory. In Conclusion we 
discuss unsolved problems and perspectives. 

2 Heisenberg double 

In this section we remind some simple results from the theory of the Heisenberg double and 
fix notations. More detailed discussion of the subject can be found in refs. |> 11- 

Let G be a matrix algebraic group and D = G x G. For definiteness we consider 
the case of the SL(N ) group. Almost all elements ( x,y ) G D can be presented in two 
equivalent forms as follows 

0 x,y) = {U J U)-\L +1 L_) = {U- 1 L + ,U- 1 L^) 

= (L+, L_)~\U, U) = ( L~ l U , LZ l U) (2.1) 

where U,U G G, the matrices L +1 L + and L_,Z_ are upper- and lower-triangular, their 
diagonal parts /+,/+ and F,F being inverse to each other: l + F = l + F = 1. 

Let all of the matrices be in the fundamental representation V of the group G (N x N 
matrices for the SL(N) group). Then the algebra of functions on the group D is generated 
by the matrix elements x l3 and y l3 . The matrices L± and U or L± and U can be considered 
as almost everywhere regular functions of x and y. Therefore, the matrix elements L± l3 
and Uij (or L± l3 and U l3 ) define another system of generators of the algebra FunD. We 
define the Poisson structure on the group D in terms of the generators L± and U as 
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follows H, I 

{U\U 2 } = 1 [r ± ,U 1 U 2 } (2.2) 


{L\,L 2 + } = 7 [r±,L\L\\ 
{Ll,L 2 _} = 1 [r ± ,L\L 2 _} 
{L\,L 2 _} = 1 [r + ,L\L 2 _] 

(2.3) 

{L\,U 2 } = 1 t + L\U 2 
{Ll,U 2 } = ir_L l _U 2 

(2.4) 

Here 7 is an arbitrary complex parameter, r± are classical r-matrices which satisfy the 
classical Yang-Baxter equation and the following relations 

[r 12 ,r 13 ] + [r 12 ,r 23 ] + [r 13 ,r 23 ] = 0 

(2.5) 

r_ = — Pr + P , r + — r_ = C 

( 2 . 6 ) 


where P is a permutation in the tensor product V (S> V (Pa ® b = b ® a) and the matrix 
C is the tensor Casimir operator of the Lie algebra of the group G. For the SL(N) group 
the solution of eqs.(2.5-2.6) looks as follows 


TV-1 


TV 


r+ = h i hi+ 2^2 e tj ® 


i =1 


i<j 


TV 


TV 


= —/ + y 

N ^ 


+ 


(2.7) 


1=1 


i<j 


&ij)ki = 5 ik 5ji and hi form an orthonomal basis of the Cartan subalgebra of the 


where 

SL(N) group: yj i(i + 1 )hi = ELi e kk - ie i+ i,*+i 
In eqs.(2.2-2.6) we use the standard notations from the theory of quantum groups 0-0 


for any matrix A acting in some space V one can construct two matrices A 1 = A® id 
and A 2 = id ® A acting in the space V ® V, and for any matrix r = Ea r i( a ) ® r 2 (a) 
acting in the space V ® V one can construct matrices r 12 = J2 a r\(a) ® r 2 (a) ® id, r l? ’ = 
E a n(a) ® id® r 2 (a) and r 23 = E a (a) ® r 2 (a) acting in the space V ® V ® V. 


The group D endowed with the Poisson structure (|2.2| - |2.4|) is called the Heisenberg 
double D\ of the group G. It is not difficult to show that the matrices L± and U have 
the same Poisson structure (2.2-2.4) and we shall need the Poisson brackets of L±, U and 
L ± , U [0 


{ Ll~L 2 } 
{LlU 2 } 
{ Li,U 2 } 
{U\U 2 } 


0 for any a,(3 = +, 
-lL l ± U 2 r± 

- 1 L\lJ 2 r ± 

0 


( 2 . 8 ) 


The cotangent bundle of the group G can be considered as a limiting case of the Heisenberg 
double . Namely, in the limit 7 —> 0 and L± —► 1 + 7 E± the Poisson structure of the 
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Heisenberg double coincides with the canonical Poisson structure of the cotangent bundle 

T*G. 

Up to now we considered the Heisenberg double of a complex Lie group. One can 
single out some real forms by means of the following (anti)-involutions: 

1. SU(N) form for imaginary 7 


U* = U~\ L\ 

2. SU(N ) form for real 7 |T3], IT 2 


3. SL(N) form for real 7 


U* = U, L* 


u * = U, L* 


LZ 1 , L*_ = L ; 1 (2.9) 

L_, L*_ = L + (2.10) 

L + , LZ = L_ (2.11) 


It can be easily checked that these (anti)-involutions are compatible with the Poisson 
structure ( 2 . 2 - 2 .4). 


3 Hamiltonian lattice Yang-Mills theory 

In this section we consider Hamiltonian lattice gauge models based on the assignment of 
the Heisenberg double to each link of the lattice. We begin with the case of an arbitrary 
graph (regular hyper-cubic lattice, triangulation of a surface, the Bruhat-Tits tree, sim- 
plicial complexes and so on) which is described by a set of vertices and a set of links. 
Each link is thought of as either a path connecting two vertices V\ and V 2 or a closed path 
with a marked vertex (tadpole). Two vertices can be connected by any finite number of 
links. Such a graph is certainly just an arbitrary connected Feynman diagram. 

Let us now consider some vicinity of a vertex v which does not contain other vertices 
and closed paths. Let us denote the paths which go from the vertex v by Ii(v),...,In v (v). 
We call such a path as a vertex path. N v is a common number of the paths and if there 
is no closed path for the vertex v then N v coincides with the number of links going from 
v to some other vertices of the lattice. With each vertex path li(v) one associates a held 
taking values in the Heisenberg double D\. This held is described by matrices U(li[y )), 
L + (li(y )) and L_(h(n)) with the Poisson structure (2.2-2.4) and helds corresponding to 
different paths have vanishing Poisson brackets. 

Let us now attach to the vertex v the following Gauss-law constraints [TJ 

G ± {v) = L±(h(v))L±(l 2 (v)) ■ • ■ L ± (l Nv (v)) = 1 (3.12) 

One can easily check that these constraints have the following Poisson brackets 

(v), G^_(u)} = 7 [r±, Gj»G 5 »] 

(Gi(u), G 2 _ (u)} = 7 [r±, Gi(u)Gl(u)] 

{Gi(u), G 2 _(v)} = 7 [r+, Gl(v)G 2 _(v)} (3.13) 

We see that these Poisson brackets vanish on the constraints surface G± = 1 and therefore 
they are hrst-class constraints. Thus one can consider gauge transformations which are 
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generated by these constraints. Namely, for any function X of the matrices U(li(v )), 
L±(k(v)) an infinitesimal gauge transformation looks as follows 

SX = {X, tr (G + (u)Mu) + G_(u)£_(u))} (3-14) 

where are gauge parameters which do not depend on U and L±. 

Using the Poisson brackets (2.4) for U and L± one can easily verify that these constraints 
generate left gauge transformations of the field U 

U{k{v)) ^ g{k{v))U{k{v)) (3.15) 

Let us note that g(lj(v)) depends not only on the vertex v but on the vertex path Uiy) 
and the fields L± as well. But the gauge transformations for different vertex paths are 
certainly not independent. 

A remarkable feature of these constraints is that they form the same quadratic Poisson 
algebra as the matrices L± do. In the limit 7 —> 0, L± —> 1 + 7 E± one gets the following 
Gauss-law constraints 

N v 

C±{v) ='£E±{l i {v)) = 0 (3.16) 

i=1 

which form the Lie algebra and were used by Kogut and Suskind |3j (more exactly they 
used the constraints C(v) = C + (v) — C_(v)). Thus the usual lattice gauge theory can be 
thought of as a particular case of the models under consideration. 

Let us remark that the choice of the constraints G±(v) is not unique. One can use 
any constraint of the form 


G±( v,a) = L ± (l n (v))L ± (l n (v)) ■ ■ -L ± (l„ (v)) = 1 


(3.17) 


where a E Symm(N v ) is some permutation of 1,2,..., 

These constraints form the same Poisson algebra ( fi. 13| ) and in the limit 7 —> 0 coincide 
with C±(v). However at finite 7 only that constraints, which differ from each other by a 
cyclic permutation, are equivalent. Thus with each vertex one can in principle associate 
(N v — 1)! nonequivalent constraints. 

Repeating the same procedure for all of the vertices one gets the phase space which is 
the direct product of the Heisenberg double over all of the vertex paths and a set of the 
Gauss-law constraints attached to the vertices. The Gauss-law constraints corresponding 
to different vertices have vanishing Poisson brackets. Taking into account that for each 
link there are two vertex paths one sees that one has placed on each link two different 
Heisenberg doubles. However one can impose on the fields attached to one link the 
following constraints 

U-\1)L ± (1) = Ll\2)U(2) (3.18) 


Comparing eq. (|3.18| ) with eq.(2.1) one concludes that the fields U( 1), L±(l) and U( 2), 
L±(2) are just different coordinates on the same Heisenberg double. Thus the phase space 
of the model is the direct product of the Heisenberg double over all links: Y[u n k s D\. 

To the moment we have described the phase space of the model and imposed the 
Gauss-law constraints . The next problem is to find gauge-invariant functions on the phase 
space and to form from them a Hamiltonian. The simplest gauge-invariant functions can 
be obtained using a well-known theorem from the theory of Poisson-Lie groups BBS 
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which states that generators of the ring of the Casimir functions of the Poisson algebra 
(2.3) have the following form 

h k = tr {L + L~_ l ) k = tr L k (3.19) 


Due to the fact that the Gauss-law constraints Q3.12| ) depend only on L± these functions 
are gauge-invariant. The functions (|3.19 ) are not the only ones gauge-invariant. Just as 
in the case of the usual lattice gauge theory one can construct Wilson line observables. 
So let us consider a loop with is formed by the oriented links Zi, Z 2 , Each oriented 

link which goes from a vertex v to a vertex u (it may be the same vertex) can be denoted 
by two vertex paths as follows 

l(v,u) = li(v)lj(u) (3.20) 


where i and j are the numbers of the vertex paths (we have numbered all vertex paths). 
Thus the loop is described by 2 n vertex paths as follows 


^1) ^2) ■■•jin hi ('D )^j2 (^2) j ^*2 (^2)^3 (^3) j ■ --i hn(,Vn)lj i (d) ( 3 - 21 ) 

i.e. we denote the link l a going from v a to v a+ i as li a (v a )lj a+1 {v a + 1 ) and G, j a = 
1) 2,..., N Va . 

In the usual lattice gauge theory the gauge-invariant Wilson line observable corresponding 
to the loop is of the form 


W(h ■ ■ • Q = tr [/(/ n (n 1 ))D(/ J 2 (n 2 )) • • ■ U(l in (v n )) (3.22) 

We are looking for a similar expression for the Wilson line observable which coincides 
with (|3.22|) in the limit 7 —> 0. It appears that to get a corresponding generalization one 
should attach to every vertex v a some vertex operator V(v a . j a , i Q ) which depends on the 
vertex paths h a (v a ) and lj a (v a ). Then the gauge-invariant Wilson line observable for a 
loop without tadpoles looks as follows 


W(h ■■■l n ) = tr U(h)V(h, i 2 )U(i 2 )V(j 3 , U(i n )V(ji,ii) (3.23) 

where we use short notations 

u(hM) = u(i a ) 

V(v a ,j a ,i a ) = V(j a ,i a ) 

The choice of the vertex operator V ( j , i ) is not unique. The simplest vertex operators 
have the form 


V(j,i) = 


L ej+ 1 {j + 1 ) ■ • • (»-l) if j < i - 1 

1 if j — i — 1 

L fj(j)--' L fi(i) if 3>i 


(3.24) 


where e k , = ±. 


To prove eq. (|3.23| ) it is enough to show that the combination 

U(l k (u))V(v,j,i)U(h(v)) = U(lj(v))V(v,j, i)U(h(v)) 


(3.25) 


is gauge-invariant under the transformations generated by the constraints G±(y). It is not 
difficult to do by using formulas (2.2-2.4, 2.8) for the Poisson brackets of U , U and L±. A 
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Hamiltonian of the model can be now written as some combination of the gauge-invariant 
functions (|3.19| ) and ( |3.23| ). 

Thus we have described the phase space, the constraints and gauge-invariant observ¬ 
ables which can be used as Hamiltonians of the models and now we are passing to a 
particular model on the regular hyper-cubic space lattice. We shall show that this model 
is the lattice-regularized Yang-Mills theory. 

So let us consider the regular hyper-cubic space lattice in d-dimensions. In this case 
an arbitrary vertex of the lattice can be denoted by a vector n = (ni, ..., n d ) with integers 
n l . We choose some orientation of the lattice and denote the orthonormal lattice vectors 
which define the orientation as e*, i — 1,..., d and introduce the notation e_; = —e*. With 
each link one can associate a positive link (n, e*) and a negative link (n + e*, — e*). The 
vertex paths of the vertex n are therefore oriented links (n, e a ), a = ±1,..., ±d. We place 
on each vertex path (n, e a ) (or oriented link) the Heisenberg double which is described 
by the fields U(n, a), L±(n, a). Due to eq. (|3.18|) we have the following relations 


U {n + e, : , -i) 
L±{ n + e i5 -i) 


= U(n,i) 
= L±(n,i) 


(3.26) 


Let us note that in the limit 7 —> 0, L± —> 1 + 7 E± one gets the usual equations 
f/(/ _1 ) = f/ _1 (/) and E(l _1 ) = —U~ l {l)E{l)U{l) for any oriented link l. 


There are (2 d — 1)! nonequivalent choices of the constraints (|3.12|) . In the paper we 
use the following Gauss-law constraints 


G±(n) = L±(n, —l)L±(n, l)L±(n, -2)L±(n, 2)---L±(n, -d)L±(n,d) 

= L±(n - ei, l)L ± (n, 1) • • ■ L ± (n - e d , d)L±(n,d) = 1 (3.27) 

We see from this expression that it is natural to introduce the notation G±(n, i) = 
L±(n,—i)L±(n,i). One can easily verify that the (anti)-involutions (2.9) and (2.11) are 
compatible with the constraints. A Hamiltonian of the model can be written in the 
following form (which is certainly not unique) 

H = Ey- d £ tr (L\t) - 1) + A1 £ («'(□) + «"(□)) (3.28) 

' links plaqettes 

where the summation is taken over all positive and negative links and over all plaqettes, e 
is the coupling constant and a is the lattice length. The Wilson term W(d) is determined 
by eqs. (|3.22|,|3.24j) and is equal to 

W(Dij) = trf/(n;z)H(n + e i ;-f,j)H(n + eqj)l/(n + e i + e j ;-j, -i) 

U{n + ei + ej; -i)V{ n + i, -j)U{ n + e^; -j)V{ n; j, i ) (3.29) 

This formula can be simplified if one uses the following equation expressing U( n + e*; —i) 
through U(n]i), L±(n,i) and L±(n,i) 

U{ n + e,; —i) = L±( n, «)f/ _1 (n, i)L±{ n, i) = L±( n + e*; —i)U~ l { n, i)L±( n, i) (3.30) 

Using this equation and eq. (|3.24|) for the vertex operators one gets for W(p) 

VU(njj) = trU(n;f)U ii (n + e i )U(n + ej; j)V^ ?: (n + ej+ ej) 

U~ 1 (n + e,-; i)^-(n + ej)L7 —x (n; j)Vji(n) (3.31) 
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where 


Vij(n) 


L ei (n,i)G ei+1 (n,i + 1) • ■■G ej _ 1 (n,j - l)L e .(n, -j) 

1 

(I'm (n, j)Gfj, j+1 (n, j + 1) • n, i - l)L w (n, -z))^ 1 


if i < j 

if i = j (3.32) 

if i > j 


Let us mention that the Wilson line observable corresponding to an arbitrary loop formed 
by oriented links li,l 2 , has the same form 


W(hl 2 ...l n ) = tr U(l l )Vi lh U(l 2 )V hh ---U(l n )Vi n i 1 (3.33) 

where one should take U(l) = f/(n;z) for any positive link l = (n, e*) and G(/ _1 ) = 
f/ _1 (n, i) for any negative link l~ l = (n + e*, —e*) and the vertex operator Vi k i k+1 is given 
by eq. (|3.32|) . Let us stress that in eq. (|3.33|) one has to use G _1 (n, i) for any negative link 
but not U( n, i). 

One can choose for example the positive sign for all of and negative sign for all of /q. 
Then the vertex operator V tJ (n) has the following transformation law with respect to the 
anti-involution (2.9) which singles out the SU(N) real form 


VS(n) = V M (n) 


(3.34) 


The formula (|3.34|) ensures that the Wilson line observable W{lil 2 ...l n ) is complex- con¬ 
jugated tO W{l n l n _ L-.il). 

Now taking into account that in the limit 7 — > 0, L± —> 1 +jE± the vertex operator Vij 
goes to the unity one recovers the Gauss-law constraints and the Hamiltonian of the usual 
lattice gauge theory ||. Thus we have shown that in the case of the regular hyper-cubic 
space lattice the model proposed is just a lattice-regularized Yang-Mills theory. 


4 Conclusion 


In this paper we considered the Hamiltonian formulation of classical gauge models on an 
arbitrary lattice. In this formulation we placed on each link the Heisenberg double of 
a Lie group and attached Gauss-law constraints to each vertex of the lattice. We have 
shown that the models on the regular hyper-cubic lattices correspond to lattice-regularized 
Yang-Mills theory. We discussed only the case of the pure Yang-Mills theory. It would be 
very interesting to include fermions to this construction. It seems to be plausible that a 
proper description of fermions would lead to a lattice version of the Faddeev-Shatashvili- 
Mickelsson 2 -cocycle [|T4|, |15[ [TO . 


Another interesting problem is to classify all integrable lattice gauge models, i.e. to 
find all corresponding graphs and Hamiltonians. 

We considered only classical theory and the next and most important problem is to 
quantize the models. There is no problem in quantizing the Poisson structure of the 
Heisenberg double. One just gets the quantized algebra of functions on the Heisenberg 
double Bing which includes as subalgebras the algebra of functions 011 the quantum 
group Fun q (G ) and the quantized universal enveloping algebra U q (Q), where Q is the Lie 
algebra of the group G. The classical r-matrices r± are to be replaced by the A-matrices 
R±(q) = 1 + ih^r± + • • •, where q = e lhl . Thus a real 7 corresponds to q lying on the 
unit circle of the complex plane and an imaginary 7 corresponds to a real q. It is of 















no problem to check that in quantum theory the Gauss-law constraints are first-class 
constraints and commute with the quantum Wilson line observables and therefore with 
the quantum Hamiltonian. Let us note that q has a nonpolynomial dependence on the 
Planck constant Ti and thus already ’’tree” correlation functions of the models will have 
a nonpolynomial dependence on h as well. It seems to be an indication that correlation 
functions of the models correspond to a summation over inhnitely-many number of the 
usual Feynman diagrams. It is not excluded that the parameter 7 plays the role of an 
infrared cut-off. Due to the fact that there is the additional parameter 7 for the models 
one may expect that these models have more rich phase structure than the usual lattice 
gauge theory. 

Let us finally notice that g-deformed lattice gauge theory was considered in refs. ||T7|, 
P~8| . p~0| . pOjl in connection with the Chern-Simons theory. 
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